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Abstract: We report on the determination of the deconfining temperature in 
SU (3) pure gauge theory, using the Symanzik tree level improved action, on 
lattices of size 3 x 12 3 , 4 x 16 3 , 5 x 20 3 , 6 x 24 3 . We find that the asymptotic 
scaling violation pattern is similar to the one observed using the Wilson ac- 
tion. We conclude that the irrelevant operators do not affect, in the range of 
couplings considered, the lattice j3 function. An analysis based on an effective 
coupling formulation shows an apparent improvement. 



1. Introduction 

In an earlier work fl|] we have reported on a high statistics investigation of 
the effect of the Symanzik tree improvement 0, based on the study of the 
asymptotic scaling properties of the deconfinement transition in the 577(2) 
pure gauge theory. In this letter we extend the analysis to the deconfinement 
transition in S77(3) theory. Our main purpose is a better understanding of the 
effect of the perturbative improvement in the context of 4 dimensional gauge 
theories. In this respect the deconfinement temperature T c is an observable 
particularly well suited to give unambiguous answers. As early noticed [§-0 
it can be accurately measured; being a bulk effect it has a slow volume depen- 
dence and most importantly it is not affected by perturbative contributions 
which violate the scaling. We determine the critical coupling of the decon- 
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finement transition in lattices with N T = 3, 4, 5, 6, and study the pattern of 
asymptotic scaling violation of T c , in comparison with the existing Wilson ac- 
tion simulations (for good reviews see |],[7) and references therein), both in the 
"bare" scheme and in the effective schemes based on the expectation values of 
plaquettes. In the bare scheme we find that the two lattice formulations give 
quite similar results, taking into account the redefinition of the A scale, which 
corresponds to the perturbative prediction whitin 10%. This is a confirmation 
that the Wilson and Symanzik actions are in the same universality class, and 
that the f3 function is not affected, in the "bare" scheme, by the irrelevant ope- 
rators connecting the two lattice formulations: on the other hand the results 
show further that with the currently accessible lattices the two-loop perturba- 
tive approximation is not sufficient to model the renormalization group flow 
in the scaling region. In sec. (fj) we give some details on our simulation and 
on the numerical results, discussing the consequences of an effective coupling 
parameterization. In sec. @ we discuss a test of our simulation, aimed to 
verify the first order nature of the deconfinement transition through a F.S.S. 
analysis on N T = 3 lattices. 



2. The simulation 



We have considered the SU(3) gauge theory on asymmetric lattices using the 
lowest order tree improved Symanzik action 

S I = ReTr (J*7 lxl - l[/ lx2 ) , (1) 

defined in order to remove O (a 2 ) scaling violations in the Callan- Symanzik 
equation for physical quantities measured on the lattice: we refer to our pre- 
ceding work m and references therein for a more detailed discussion of the 
improvement strategies and of the results existing in literature. We have per- 
formed simulations on lattices N T X N% = 3 X 12 3 , 4 x 16 3 , 5 x 20 3 , 6 x 24 3 ; 
the update is done with a combination of Heat-Bath and of overrelaxation 
sweeps applied to the SU(2) subgroups ||. All the data have been obtained 
on the APE machines installed in Pisa, including a 128 processors machine 
with a peak speed of 6 GigaFlops. We have measured the order parameter of 
the deconfining transition, the spatial average of the Polyakov line 

N t 

P = N^Y, T *U U *& t ) ■ ( 2 ) 
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As usual (see for example H) to consider the projection of this quantity on 
the nearest Z (3) axis, 



ReP argP G (-§ , f 

RePexp(-2ivr/3) argP e (f, vr) . (3) 
ReP exp (2 % tt/3) argP e ( — 7r, 



in order to avoid the vanishing of the order parameter, on finite lattices, in- 
duced by the tunnelling between broken phases. We locate the deconfinement 
transition by studying the Polyakov susceptibility 

x(n) = Jv*((n) a -<n> 2 ) (4) 

which exhibits a clear peak structure, and the specific heat of the lattice action. 
The density of states method [|Hj] has been used to extrapolate data coming 
from runs at different (3 values and to combine them in a single curve: the data 
autocorrelation has been taken in account using a standard binning procedure, 
varying the bin size and testing the existence of a plateau in the variance of 
the binned quantities. In order to reduce the bias induced by the finite sample 



size we have used the techniques discussed in ||11|| : "standard" and "jackknife" 

estimators have been defined, combining them with coefficients appropriate to 

cancel the term in the bias. In fig. ([!]) is given an example of the resumma- Fig.|1| 

tion of the Polyakov susceptibility around the deconfining transition on the 

N T — 3 lattice, where the actual simulation values are addressed by a black 

square. A more extended discussion of the numerical simulation will be given 

elsewhere. In tab. ([I]) we report the values of the (pseudo) critical coupling 

on the different lattices, obtained by fitting the peaks with a gaussian curve, 

and the corresponding critical temperature in units of Aj, obtained by means 

of the two loop R.G. formula 



T r 1 




(5) 
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They are also reported the corresponding quantities obtained in simulations 



using the Wilson action |3]-|^|12| . In the last column the corresponding ratio of 
the A parameters in the improved and standard formulation is listed, showing 
at N T = 6 a 10% discrepancy from the perturbative prediction 

A/ 5.29210(1) . (6) 
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lt is instructive to compare different actions and schemes by normalizing 
the results in order to have the values at smaller a coincide, as in fig. (0). Fig. [2] 
It is quite apparent that the asymptotic scaling violations for the Symanzik 
and Wilson action, in the "bare" scheme, are similar, as already noted in the 
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Fig. 1. Susceptibility of the Polyakov line on a 3 x 12 lattice 



Table 1 

Pc values on different lattices 



N T 


Pc (Improved) 


r c /Aj 


Pc (Wilson) 


T c /A w 




3 


3.90812(7) 


13.937(2) 


5.55(1) 


85.70(96) 


6.39(7) 


4 


4.07252(13) 


12.505(2) 


5.6925(2) 


75.41(2) 


6.030(2) 


5 


4.19963(14) 


11.497(2) 


5.7933(3) 


68.5(1) 


5.95(1) 


6 


4.31466(24) 


10.870(3) 


5.8941(5) 


63.05(4) 


5.800(5) 



SU(2) simulation [|]]. In the same figure we show the same data parameterized 
in terms of the expectation values of plaquettes, in the manner of Parisi ||14|| : 
to define effective couplings, the perturbative predictions on the lxl and 
1x2 plaquettes [O] have been used, together with the results of simulations 
on large symmetric lattices. As already noted @ for the results obtained with 
the Wilson action, this parameterization shows an apparent asymptotic scaling 
already for N T > 4; our simulation with the Symanzik tree improved action 
shows an apparent scaling already for N T — 3. 
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Fig. 2. Asymptotic scaling figure for the critical temperature 

We recall that the Symanzik approach is based on perturbation theory: hence 
the fact that the improvement is apparent when an effective coupling is used 
seems a confirmation of the work of Lepage and Mackenzie []I5|| , which asserts 
that the bare coupling is a bad expansion parameter, and that the perturbative 
series, when properly reorganized, work well already on the accessible lattices. 
In other words the effect of the Symanzik improvement is "buried" by the 
large renormalizations induced by the lattice formulation: the improvement 
shows up when a scheme is chosen which minimizes these finite corrections. 



3. Finite size study for N T = 3 



In order to control at least in part the systematic error connected to the 
finite size of the lattices, we have performed for iV T = 3 an additional set of 
simulations with N a = 6, 7, 8, 9, 10, 11, to test the dependence of the critical 
coupling on the volume of the system, and to check the expected scaling 
relations for a first order phase transition. We report in tab. (||D the values of 
the position of the peak of the susceptibility of the Polyakov line and its height, 
resulting from a gaussian fit with a self consistent data windowing around the 
peak. In fig. (^) we compare the results for the Polyakov line susceptibility: we Fig. |3| 
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Fig. 3. Polyakov susceptibility for N T = 3, N a = 6, ... 12 



Table 2 



Results of the finite 


size study on the N T = 3 lattices. 




^max 




6 


3.90588(34) 


0.00568(5) 


7 


3.90525(53) 


0.00514(12) 


8 


3.90825(46) 


0.00505(12) 


9 


3.90859(24) 


0.00450(12) 


10 


3.90765(5) 


0.00415(16) 


11 


3.90853(8) 


0.00436(9) 


12 


3.90812(7) 


0.00416(7) 



note that the scaling of the peak is compatible with the characteristic behavior 
in presence of a first order phase transition 

x (n) * K . (7) 
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Most importantly for our purposes, we note that the shift of the pseudo- 
critical coupling corresponding to the peak of the transition is modest, giving 
confidence in the values and the errors quoted in tab. ([!]) for the lattices with 
iV r > 3. 



4. Conclusions 

In this work we have performed a study of the deconfinement phase transition 
in the pure gauge £77(3) theory, using the Symanzik tree improved action and 
determining the critical couplings on lattices with N T = 3, 4, 5, 6. The main 
conclusion is that, as in the case of the deconfinement transition in SU(2) 
theory Jl| , the same pattern of asymptotic scaling violation is observed. The 
asymptotic scaling violations are much reduced if the two loop perturbative 
expression is used in conjunction with effective couplings, based on UV domi- 
nated quantities like the plaquette expectation values, and comparison with 
the results of simulations with the Wilson action, in this "effective" schemes, 
shows a slighlty precocious scaling, which is tempting to interpret as an effect 
of the improvement. Although we do not have under control the systemati- 
cal errors due to the finite volume, the analysis of the finite size behavior in 
the N T = 3 lattices makes us confident that the quoted results are not much 
affected by the finite size shift of the transition. 
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